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Abstract. Ions stored in Penning traps may have useful applications in the field 
of quantum information processing. There are, however, difficulties associated 
with the laser cooling of one of the radial motions of ions in these traps, namely 
the magnetron motion. The application of a small radio-frequency quadrupolar 
electric potential resonant with the sum of the two radial motional frequencies has 
been shown to couple these motions and to lead to more efficient laser cooling. We 
present an analytical model that enables us to determine laser cooling rates in the 
presence of such an 'axializing' field. It is found that this field leads to an averaging 
of the laser cooling rates for the two motions and hence improves the overall laser 
cooling efficiency. The model also predicts shifts in the motional frequencies due 
to the axializing field that are in qualitative agreement with those measured in 
recent experiments. It is possible to determine laser cooling rates experimentally 
by studying the phase response of the cooled ions to a near resonant excitation 
field. Using the model developed in this paper, we study the expected phase 
response when an axializing field is present. 
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1. Introduction 

Quantum information processing (QIP) has been the subject of much research in 
recent years because of its potential ability to provide a polynomial or even exponential 
speedup over conventional computing when solving certain problems, perhaps most 
notably the factoring of large integers, the searching of unsorted lists and the 
simulation of quantum systems [1, 2]. A key requirement for any system that is 
to form the basis of a quantum computer is that it must be well isolated from its 
environment in order to prevent decoherence. For this and other reasons, ions trapped 
in radiofrcquency (RF) traps have proved to be an excellent system and many of the 
basic requirements for QIP have been demonstrated using them [3, 4, 5, 6, 7, 8]. 
However, the large amplitude oscillating fields used for confinement in RF traps mean 
that ion heating rates can be rather large and can lead to a loss of coherence during 
gate operations [9, 10]. With its purely static electric and magnetic fields, the Penning 
trap may therefore provide an attractive alternative system for research towards QIP. 
The operation of a universal quantum computer based on quantum gates is extremely 
challenging since it requires very high fidelity and low decoherence. On the other hand, 
direct quantum simulation (DQS) in which trapped ions are used to directly simulate 
other quantum systems may prove to be a more tractable experimental prospect. An 
array of ions in a Penning trap has recently been proposed as an ideal system in 
which to perform DQS of a range of systems of interest to the condensed matter 
community [11, 12, 13]. Work is also being carried out on the use of novel Penning 
trap designs for quantum information processing [14, 15]. 

By using a Penning trap located in a superconducting magnet it is also possible 
to reduce the influence of fluctuations in the ambient magnetic field, which currently 
limit the coherence times that are observed in many RF trap experiments. In order to 
combat the effects of fluctuating magnetic fields there has been much recent interest 
in the use of first-order magnetic field free transitions as qubit transitions in QIP 
studies. These are transitions between two internal electronic states whose Zeeman 
shift gradients are at some point equal, so that the transition frequency is (to first 
order) insensitive to changes in the magnetic field [7] . In many cases these transitions 
occur at rather large magnetic fields, making them well suited for use with ions trapped 
in a Penning trap [16]. 

The two key difficulties with working in the Penning trap are associated with 
the laser cooling of the ions. Firstly, the large magnetic field used for trapping 
(typically a few tesla) leads to rather large Zeeman splittings and so the number 
of transitions that must be addressed to form a closed cooling cycle is greater than 
in a RF trap. This increases the technical complexity of an experiment, but has not 
proved insurmountable [17, 18]. The second difficulty is that the ions in a Penning trap 
have a rather more complex motion that includes an unstable component, called the 
magnetron motion. Under normal circumstances laser cooling of this component of the 
motion is necessarily inefficient. A technique called axialization has been developed 
and demonstrated to improve the efficiency of laser cooling in the Penning trap by 
coupling the two radial motions [19, 20]. An analytical model for axialization is 
presented in this paper. This model exhibits a number of features that have been 
observed in recent experimental work on the axialization technique. 

We first review both ion motion in a Penning trap and a simple analytical model 
for laser cooling of the two radial motions. We then proceed to introduce the effect 
of the coupling of these two motions into the model and derive damping rates in the 
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Figure 1. (a) Electrode structure of an idealised Penning trap, (b) Non- ideal 
Penning trap with a split ring electrode. This is similar to a trap used in recent 
experiments. Note that a section from the left image and one of the ring segments 
from the right image have been removed for clarity. 



presence of axialization. Finally, we examine how an external excitation of one of the 
motions can reveal information on the damping rates and the coupling rate between 
the motions. 

2. Motion in the Penning trap 

A detailed description of the Penning trap can be found in, for example, [21, 22]. 
In brief, an idealised Penning trap consists of electrodes whose inner surfaces are 
hyperbolae of revolution about the z axis (figure 1(a)). For positively charged ions 
the two endcaps are given a positive bias with respect to the ring electrode. This 
leads to a harmonic electric potential <p that confines in the z direction, but results 
in an outward force in the radial plane. Radial confinement is achieved by applying 
an additional static magnetic field of typically a few tesla along the axis of the trap. 
This forces the ions into cyclotron-like loops. In practice it is often convenient to use 
electrode shapes that differ somewhat from the ideal. Furthermore, in order to apply 
the axialization technique described in this paper it is necessary for the ring electrode 
to be split into four radial segments (see figure 1(b)). This enables an additional small 
radially quadrupolar electric potential to be applied. 

The force on an ion with charge e due to the trapping field is given by the Lorentz 
force 



where B = —Bz and <f> = Uq (2z 2 — x 2 — y 2 ) / (2z 2 + rjj). In the z direction this 
results in simple harmonic motion with a frequency 



F = ev x B - eV<^, 



(1) 




(2) 
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In the x and y directions (1) leads to the coupled equations 
x + LU c y —x = 0, 

LU 2 

V - v c x - -^y = 0. (3) 

Here oj c is the true cyclotron frequency eB/m. These equations can be solved by 
defining a new variable u — x + iy, so that 

uj 2 

u — iuj c u ^-u = 0. (4) 

Trying a solution corresponding to circular motion, u = uoe lwt , leads to 

u e wt (uJ 2 - Wc td + = 0. (5) 



2 



This has roots 



2 V 4 2 

LOr 

= y±<"i> (6) 

where u>\ has been defined to equal the square root term. Thus the radial motion 
is a superposition of two circular motions, the rapid modified cyclotron motion with 
frequency lo' c = lu c /2 + u>i and the usually much slower magnetron motion with a 
frequency u> m = u> c /2 — u)±. 



3. Laser cooling in the Penning trap 

The model presented in this paper is an extension of work first presented in [23]. Only 
a brief description of the key elements is given here. The model deals with the case 
of a single ion being Doppler cooled by a laser with a frequency z^l, detuned from 
some strong transition at a frequency z/ . The laser beam is directed parallel to the 
x axis, but can be offset by an amount yg in the y direction (figure 2). Absorption 
of photons from the laser beam leads to a net force on the ion in the x direction. 
Photon re-emission is assumed to occur in random directions and hence to cancel to 
zero. Since the re-emission process is important in determining the ion temperatures 
achievable through laser cooling, it follows that the model as it stands cannot be used 
to calculate limiting temperatures. Instead it yields the rates at which ion motions 
are laser cooled. 

The Gaussian intensity profile of the beam means that there is a variation in the 
photon scattering rate, and therefore the force on the ion due to the laser, dependent 
on the position of the ion in the y direction. Similarly, the velocity of the ion in the x 
direction leads to a Doppler shift of the laser frequency and hence the ion experiences 
different regions of the Lorentzian lineshape of the cooling transition as it moves in 
the trap. The force on the ion is therefore dependent also on x. Providing the ion is 
cold, so that the amplitude of its motion and its velocity are small, the variation in 
photon scattering rate with both y and x is approximately linear (figure 3). 

The total force on the ion due to its interaction with the laser beam is in the x 
direction and can be written as F = Fq + F a + Fp. Here Fq is a constant that leads 
only to a slight offset in the equilibrium position of the ion and is henceforth ignored. 
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Figure 2. Diagram showing the experimental situation assumed in the laser 
cooling model. The Gaussian on the left shows the intensity of the beam as a 
function of y and indicates the changing intensity experienced by the ion as it 
moves counter-clockwise in the trap. 




Figure 3. Left: If the amplitude of an ion's motion is small compared to the waist 
of the beam, the variation in scattering rate R as a function of y is approximately 
linear. Right: Similarly, if the ion velocity is small then the Doppler shift as 
the ion oscillates is small compared to the linewidth of the transition and so the 
change in scattering rate as a function of x is approximately linear. 



F a and Fp are the components of the force that are considered to vary linearly with 
ion position y and velocity x respectively The constants of proportionality a and (3 
arc defined such that F a = —2amy and Fp = —2[3mx. The values of a and /3 are 
determined by calculating the slopes of the power-broadened scattering rate curves 
shown in figure 3 at the points where y and x are zero. Note that a and (3 can be 
changed by varying the beam offset, yo, and frequency, z/l, but the two parameters 
are not independent. 

Adding these forces to the equations of motion obtained in section 2 yields the 
equations of motion in the presence of laser cooling 

x + uj c y + 2(3x - -Z-x + 2ay = 0, 

y-u) c x- -^-y = 0. (7) 
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Assuming that (3 <C uj c and a <C uj 2 /2, corresponding physically to the assumption 
that the laser interaction is sufficiently weak that it does not significantly affect the 
ion oscillation frequencies, and changing as before to a complex variable u = x + iy, 
leads to 

u + {(3 - iu c ) u - + iaj u = 0. (8) 

Trying once again solutions corresponding to circular motion of the form u = uoe lwt , 
where uj is now also complex, gives 



u> 2 



u! 2 — (u> c + i/3) uj + -2- + ia 



= 0. (9) 



This has the roots 



oj = — - — ±du>f- — + — ia. (10) 

A Taylor expansion of the square root term with the assumption that a <IC oof and 
(3 <C lo\ leads to the solutions 

, . (3uj' c -a 

and 

CJ=LJ m + l — . (12) 

ZU>\ 

The frequencies of motion, given by the real part of uj, are therefore the same as those 
obtained in the absence of laser cooling. 

Looking at the damping rates, given by the imaginary parts of the solutions, 
the difficulty in achieving efficient laser cooling in the Penning trap is immediately 
apparent. Positive damping of the modified cyclotron motion requires that (3ui' c > a, 
whereas cooling of the magnetron motion requires that a > /3u) m . It should be 
remembered that a and f3 are related by the overall scattering rate — offsetting 
the beam further may increase the magnitude of a, but will always reduce the overall 
scattering rate and therefore (3. This means that to obtain efficient cooling of the 
modified cyclotron motion a red-detuned laser (positive (3) is needed, with perhaps 
a slight offset of the beam in the y direction (negative a). Strong cooling of the 
magnetron motion requires a blue-detuned laser (negative (3) with a small offset in 
the — y direction. These requirements are incompatible and therefore it is impossible 
to achieve efficient laser cooling of both motions simultaneously in a conventional 
Penning trap. Limited damping of both motions can, however, be achieved if the 
criterion uj m < a/ (3 < uj' c is satisfied, where both a and j3 must be positive [24, 25]. 
This corresponds to a red-detuned laser and an offset of the beam such that it 
propagates in the same direction as the ions travel around the centre of the trap. This 
is confirmed by numerical calculations of the cooling rates for the modified cyclotron 
and magnetron motions as a function of the beam offset j/o and detuning z/l — ^o, and 
by numerous experiments [26, 27]. 

Figure 4 shows the results of such a calculation for a calcium ion with trapping 
parameters that are typical of those in recent experiments [18]. A beam waist of 50/im 
is assumed and the peak beam intensity used is 0.1 times the saturation intensity of 
the transition. The only way to achieve cooling of both motions is to red-detune and 
offset the laser beam, as indicated in figure 4. 
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i/ L - v I MHz z/l - i/ I MHz 

Figure 4. Modified cyclotron (left) and magnetron (right) motion cooling rates 
calculated using the described cooling model. A beam waist of 50 fim and a peak 
intensity of 0.1 X J sa t is assumed. Lighter shades indicate larger cooling rates and 
in each case the zero cooling rate contour is that which bisects the plot. The small 
ring indicates a position at which some cooling of both motions can be achieved. 



4. A model for axialization 

The difficulty in efficiently laser cooling the magnetron motion can be eased by 
applying a weak radially quadrupolar electric potential that oscillates at the true 
cyclotron frequency, lo c . This is the sum of the two radial motional frequencies and 
hence the potential acts to couple these two components of the motion together. 
The relative amplitudes of the two components are then linked so that the efficient 
cooling of the modified cyclotron motion also leads to a reduction in amplitude of the 
magnetron motion. Note that the sum frequency, rather than the difference, is used 
because the magnetron motion is unstable — energy must be added to this motion in 
order to reduce its amplitude. This is due to the fact that the saddle point of (j> in (f ) 
has a maximum in the radial plane at the origin. Throughout this paper we adopt 
a pragmatic definition of 'cooling' of the magnetron motion as being a process that 
reduces its amplitude. 

Axialization, or 'radial centering', has been used in combination with resistive 
cooling [28] and buffer gas cooling using both neutral and charged species [29, 30]. 
More recently it has been used with laser cooling [f9, 20]. Outside of the laser cooling 
community the axialization technique is also known as 'sideband cooling'. This is not 
related to the laser cooling technique of the same name [31]. Theoretical descriptions 
of axialization have been given in [32, 33, 34, 35, 36]. The analytical model for laser 
cooling in the Penning trap presented in the previous section can be extended to 
include the effect of an axializing field, and doing so helps explain various effects 
observed in recent experiments. 

Consider a radially quadrupolar potential oscillating with a frequency w a , 

? 9 

x — y 

</>a = Vq 2 COS (uj a t) . (13) 

Zr 
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In polar coordinates this becomes 

Vo r 2 , o 

^ = ( cos 9 

z i 



9) cos (ujJ) 



Vo r 2 

-f— (cos (29 + uj a t) + cos (26» - uj a t)) . 
4 r 



(14) 



which is equivalent to two counter-rotating quadrupolar potentials. Transforming into 
a frame rotating at a frequency w r by making the substitutions 9' = 9 — uj r t and r' = r, 
means that in this frame 

Vn r' 2 

<j>' a = -f— (cos (29' + (2uj r + uj a ) t) + cos (29' + (2lu, - w a ) t)) . (15) 
4 r 

Clearly if uj t = w a /2 then one of these components is static and the other rotates 
with frequency — 4uj r . It is anticipated that the oscillating quadrupolar potential will 
give rise to a coupling between motions only when close to resonance with a sum 
or difference of motional frequencies. In a frame rotating at lu c /2 the two motional 
frequencies arc ±a->i, so the sum frequency is zero. Thus the static component of 
a will be resonant, but the (rapidly) rotating component will in general not be. It 
is therefore possible to ignore the rotating component and calculate the acceleration 
in the rotating frame due to the static component alone. If x' and y' are cartesian 
coordinates in the frame rotating with a frequency uj r = u a /2 then 

^ = ^' 2 cos (29') 



v ' 



(16) 



The additional forces on a positively charged ion can then be calculated as 

eV , 



iv = - 



2r 2 ' 
z 'o 



F y ,= 



eV ( 
2r 



Making the substitution v 



F v = R, 



i 


x'- 



iF v 



iy' allows this to be written as 
2r 2 V • 



(17) 



(18) 



By introducing e = eVg/2mrQ, the force due to axialization is simply given by —emu*. 
This force can now be included in the laser cooling model. 

Equation (8), which describes the acceleration due to laser cooling in the u 
(laboratory) frame, is first transformed into the v (rotating) frame in which the 
axializing force has been determined. Using the relation u — vo luJlt leads to 



v + ((3 — iuj c + 2itu r ) v 



+ ia + cJ r — uj t (uj c + i(3) I v 



Adding in the axialization term then gives 
v + (/? — iuj c + 2iw-[) v 



-2- + ia + lj 2 — w x (w c + if3) ) v + ev* = 0. 



(19) 



(20) 
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The appropriate trial solution now corresponds to elliptical motion with uniform 
decay, because the rotational symmetry is broken by the v* term. This elliptical 
motion is a superposition of two circular motions in opposite senses with different 
amplitudes. In general. 

v = Ae iwt + Be- i "* t , (21) 

For each solution u>, A and B are the counterclockwise and clockwise components of 
a single motion. Note the use of lo* in the second term — the damping rate (real part 
of the exponent) must be the same for both components even though they rotate in 
different senses. 

A complete solution to (20) using this trial solution can be found in Appendix A. 
In it we define A to be half the detuning of the axialization drive from the cyclotron 
frequency, so that lo & = lo c + 2A. Thus the rotation frequency of the frame in which 
the axializing quadrupolar drive is static is lo t = lu c /2 + A. We express the solutions lo 
in terms of the unperturbed motional frequencies lo\ and additional (real) frequency 
shift and damping terms So and 70 such that lo = u>\ + So + i-fo- The solutions obtained 
are stated here as 



and 



where 



and 



5 Q = ±^ F \/N+^W+^^ (22) 
V2 

P AM 

M = (2a - (3lo c ) /2lo 1 (24) 

Note that the parameter M is a measure of the strength of the laser cooling. N, on 
the other hand, is dependent not only on the laser parameters, but also the amplitude 
and detuning of the axialization drive. 

Recalling that the motional frequencies in the laboratory frame are given by 
u = Lu r ± (u>i + So) and that ui T = tu c /2 + A it is apparent that having two solutions 
for Sq leads to components in the motion at four frequencies, given by 

w = uj' c + A + 5 (26) 

and 

lo = Lj m + A - So (27) 
where So takes the two values given by (22). 



5. Results of the model 



The behaviour predicted by the model depends on the relative values of the terms that 
make up the parameter TV, defined by (25). Clearly if the amplitude of the axializing 
quadrupolar drive is set to zero then the model should yield unshifted frequencies in 
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the laboratory frame, i.e. w' c and uj m . This is confirmed by setting e to zero. In this 
case N = A 2 - M 2 /4 and (22) reduces to 

S = ±A. (28) 

In the laboratory frame the frequency shifts for each motion (given by So + A) are 
therefore +2 A or zero. As we show in section 5.4, the amplitudes of the former solution 
turn out to be zero, leaving us, as expected, with unshifted frequencies. We can obtain 
the damping rates by putting 6q = ±A into (23), giving 

/3±M 



7o 



8lu c 



2ui 1 



(29) 



so either 



or 



70 = (30) 

7o = — -z • (31) 

2wi 

For the modified cyclotron motion it is the latter solution that corresponds to the 
unshifted motional frequency in the laboratory frame. The — sign in (27) means that 
the solutions for <5o are reversed for the magnetron motion, and hence it is the former 
damping rate that applies to the unaltered frequency. Thus the damping rate obtained 
for each motion is the same as that from the model without axialization. 

Interestingly, when there is some non-zero amplitude for the axialising quadru- 
polar drive, e, the behaviour of the system in terms of frequency shifts becomes 
dependent on the laser parameters. The critical factor is the balance of the three 
terms in the expression for N . Note that if the detuning of the axializing drive is 
made sufficiently large the A 2 term in N will always dominate and So will reduce 
once again to ±A, as it did for the case of no drive at all. This demonstrates that 
axialization is very much a resonant effect and the drive frequency must be close to 
the true cyclotron frequency if it is to have an effect. The effective width of the 
resonance is related to the amplitude of the quadrupolar drive and the strength of the 
laser cooling. Close to resonance A is small and the frequency behaviour depends on 
the relative magnitudes of e/u]\ and M. We shall now consider in turn three cases 
corresponding to (e/uji) 2 < M 2 , {e/wif ~ M 2 and {t/uJif > M 2 . 



5.1. Case 1: (e/wi) 2 < M 2 

In this regime N ~ A 2 — M 2 /4 and the results are to good approximation the same 
as in the previously considered case of no axialization. Calculated motional frequency 
shifts in the laboratory frame are plotted for the magnetron motion in figure 5(a), 
along with the corresponding damping (i.e. laser cooling) rates. The motion which has 
nearly unshifted frequency has the damping characteristics of the magnetron motion 
in an unperturbed Penning trap. The other motion is the modified cyclotron motion 
'dressed' by the axializing field. It has the large damping characteristic of the modified 
cyclotron motion but, as we shall see in section 5.4, we expect its amplitude in this 
regime to be very small. There is no visible coupling of the two damping rates. Note 



Laser cooling in the Penning trap 
(a) 




.Q -0.5o' 1 1 1 ' 1 1 1 1-0.2 

TO -0.50 -0.25 0.00 0.25 0.50 

lab. frame axialization drive detuning, 2A / kHz 



(c) 




CD -0.50 -0.25 0.00 0.25 0.50 

lab. frame axialization drive detuning, 2A / kHz 



11 



(1)) 




lab. frame axialization drive detuning, 2A / kHz 



(d) 




lab. frame axialization drive detuning, 2A / kHz 



Figure 5. Motional frequency shifts (solid lines) and cooling rates (dashed lines) 
relative to the magnetron motion in the laboratory frame for the cases where (a) 
(e/wi) 2 « M 2 , (b) (e/a;i) 2 = A/ 2 , (c) (e/wi) 2 > M 2 and (d) (e/a>i) 2 > A'/ 2 . 
The equivalent situation for the modified cyclotron motion differs only in that 
the cooling rate curves are swapped around. Thin solid lines have been used 
to indicate when one of the components of the motion has less than half the 
amplitude of the other, as determined from an analysis described in section 5.4. 
Plotted parameters are for Ca+ at 0.98 tesla, such that ui c = 380 X 2n kHz and 
U ) 1 = 165 X 2?r kHz. In each case |M| 2 = 0.01 kHz and a//3 = 100 X 2tt kHz. The 
inset plots are enlargements of the central region in which the behaviour for small 
axialization drive detuning can be more clearly seen. Note the change in scale in 
plot (d). 



the swapping of the curves at A = 0. The differently shaded curves correspond to 
taking either the positive or negative solution for Sq and there is a sudden change 
in which root represents which motion. Note also that a similar plot would be seen 
near the frequency of the modified cyclotron motion, in which an equivalent additional 
frequency would be seen corresponding to the 'dressed' magnetron motion. 

5.2. Case 2: (e/wi) 2 - M 2 

In this regime more interesting effects begin to occur. Although the A 2 term always 
dominates the expression for N (except for very small A), it is not always the dominant 
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term in the equation for <5 - If N = A then 



1 



S = ± — A 2 + VA 2 (A 2 + M 2 ). (32) 
In the region where A is small compared to M this expression becomes 



± -^=^A (A + VA 2 + M 2 ) 

^ (33) 



Therefore for small A the shift in motional frequencies will vary as ±v A in the rotating 
frame. As A becomes larger compared with M the curves will again tend to ±A. 

Figure 5(b) is a plot of the frequency shifts and cooling rates calculated with 
(e/u>i) = M 2 . With the exception of e, all the parameters are the same as those used 
to plot figure 5(a). The yA behaviour near resonance that was predicted to occur 
in the rotating frame is modified by the transformation into the laboratory frame. 
Further from resonance the frequencies tend to those observed in the low axialization 
amplitude regime described in section 5.1. Note that whenever (e/oji) 2 < M 2 the 
frequency shift Sq is equal to zero when the axialization drive detuning is zero. 

Looking at the damping rate curves, the useful property of the axialization 
technique is seen for the first time. If {e/ioi} 2 is nearly equal to or greater than 
M 2 , then near resonance the damping rates of the two motions converge and there is 
an intermediate damping rate for both motions. The benefit of coupling together the 
magnetron and modified cyclotron motions like this is not just that the magnetron 
damping is increased at the expense of the modified cyclotron damping, but that it is 
no longer necessary to maintain a large beam offset parameter a and small frequency 
gradient parameter (3 in order to prevent heating of the magnetron motion. The 
average cooling rate of the two motions is 



1 (0 AM /3 AM\ 
7avc " 2 U ± ~25o~ + 2 T "2*T ) 

= £ 
2' 



(34) 



The a dependence has gone and the optimum case is represented by maximising (3. 
Furthermore, reducing the beam offset increases the beam intensity seen by the ions 
and so increases the overall scattering rate and hence j3, yielding still higher mean 
cooling rates. The only disadvantage to this is that it is possible to realise a situation 
in which one of the cooling rates is negative in the absence of axialization, and so the 
frequency of the axialization drive becomes even more important. If the laser cooling 
is set to give large average damping rates but the axialization drive is not properly 
tuned to resonance, then the size of the magnetron orbit will be rapidly increased. 
This is highlighted in figure 6, which shows the cooling rate curves for the situation 
in which (3 is slightly greater than in figure 5. The ratio a/ (3 is now 10 x 2tt kHz so 
that the cooling criterion in the absence of axialization, uj m < a//3 < u' c , is no longer 
satisfied. Clearly the damping rates when the axialization drive is near resonance are 
larger, but if the drive is detuned by more than about 0.7 kHz then the damping rate 
of one component of the motion becomes negative. 
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Figure 6. Cooling rates for the magnetron motion when the ratio a/0 is chosen 
to be 10 X 2n kHz — ten times smaller than in figure 5. Although the mean 
cooling rate is larger near resonance, if the axialization drive is detuned by more 
than about 0.7 kHz one of the cooling rates becomes negative. As in figure 5(b), 
(e/uii) 2 = M 2 with both terms now being 0.1 kHz 2 . As in the previous figures, 
uj c = 380 X 2tc kHz and u\ = 165 X 2tt kHz. 



For the special case where {e/uJi) is exactly equal to Af 2 , the shapes of the 
cooling rate curves near resonance are given by substituting So = ±\/\AM/2\ into the 
equation for 70. This yields 

5.3. Case 3: (e/wi) 2 > M 2 

In this regime there is a similar damping rate coupling to that seen in the previous 
case, although the range of detunings in which the effect is significant is larger. What is 
different, however, is that as soon as (e/wi) 2 > M 2 the frequency shifts no longer cross 
at zero. Instead there is an avoided crossing — the frequency shifts approach each 
other as before for large A, but near resonance they split apart. The damping rates 
still swap over at A = as required, though at no point do the motions share the same 
frequency. If (e/uji) 2 is only slightly greater than M 2 then there is an abrupt change 
in the gradients of the frequency shifts near the avoided crossing. An example of this 
is given in figure 5(c), which shows the frequency shifts and cooling rates calculated 
using the same parameters as in the previous cases, but with (e/uji) 2 = 1.05A/ 2 . As 
(e/lui) 2 is increased further the frequencies of the two components of the motion split 
further apart and the gradients of the frequency shifts vary more smoothly. This case 
is illustrated in figure 5(d), in which the frequency shifts and cooling rates are shown 
for the case where (e/cJi) = 100M 2 . 

If (e/uJi) 2 M 2 then the approximation N = A 2 + e 2 / (4w 2 ) can be made. The 
frequency shifts in the rotating frame are then given by 



AM 



(35) 
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= ±Vn 

= ± ^ + 4 (36) 

The frequency shifts at A = are simply equal to e/ (2wi) and the hyperbolic form 
of this function is in excellent agreement with recent experimental measurements. 

The situation discussed here, where the laser cooling forces are weak compared 
to the forces due to the axialization drive, has been considered before using a different 
approach [19, 20]. The physical interpretation is that there is a continual exchange of 
amplitude back and forth between the two motions. Thus if the amplitude of one of 
the motions can be reduced by laser cooling then with axialization the other motion 
will also be damped. 

Avoided crossings resulting from the coupling together of motional modes of ions 
have been reported before by Cornell et al. [37] and Haffncr et al. [28] using a Fourier 
Transform Ion Cyclotron Resonance technique. Both groups coupled a radial motion 
to the axial motion in an attempt to extend resistive cooling of ions to the radial 
motions. By a different method Cornell et al. derived the same equation formulated 
here for the frequency shifts in the regime where the coupling rate is stronger than 
the cooling rate. The result here extends theirs into the region of weaker coupling. 
Furthermore the work by Cornell et al. did not consider the phase behaviour of the 
ions' motion. This motional phase behaviour forms an integral part of experimental 
techniques often used to measure ion cooling rates. Since such measurements make 
it possible to determine the efficacy of axialization experimentally, we shall now 
investigate the motional phase behaviour in this model. We shall also look at the 
relative amplitudes of the different motions by considering their ability to couple to 
an external driving field. 



5.4- Amplitude and phase response to a driving field 

Thus far the motional frequencies and damping rates for an ion being laser cooled in 
the Penning trap in the presence of axialization have been calculated. It was found, 
however, that the extra motional frequencies arising from the axializing field exist as 
solutions even if the amplitude of the drive, e, is zero and it was suggested that in this 
particular case these components of the motions would always have zero amplitude. 
The solutions obtained (in the rotating frame) were of the form 

v = Ae iut + Be- W ' 1 . (37) 

The real and imaginary parts of to correspond to the motional frequency and damping 
rate respectively. The relative amplitudes of the two components of each motion can be 
determined by finding an expression for the ratio of A to B. It is, however, interesting 
to derive a more general result that also includes the effect of an additional external 
force on the ion, since such a force is used in an experimental technique for measuring 
laser cooling rates [38] . 

We consider an additional dipolar field oscillating at a frequency close to that 
of one of the characteristic ion motions in the laboratory frame, so that the ion is 
driven at this frequency. The amplitude of the motion depends on how closely the 
driving frequency is tuned to one of the ion's natural motional frequencies. When 
the drive frequency is exactly resonant with a natural frequency of the ion's motion 
the amplitude of the response can be extremely large. The phase of the ion's motion 
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relative to the driving field also depends on the drive frequency — the ion will change 
smoothly from moving in phase with the drive to miving in antiphase with it as the 
drive is tuned from below to above resonance. The width of this resonance in the 
response of an ion to a driving field is directly related to the damping present in the 
system. Measurement of this width therefore allows the experimental determination 
of laser cooling rates [38] . 

A complete solution to the equations of motion of a driven, laser-cooled ion in 
the presence of axialization is given in Appendix B. It adds to (20), the equation of 
motion for the ion, an additional driving force term Fe luJt . Using the trial solution 
given in (37), where to is now the frequency of the excitation drive in the rotating 
frame, the motional amplitudes A and B can be determined as 

A= t-Q—s, (38) 

C m C c -|ef 

Fe , s 

B = —5 , 39 

where 

C m = -26wi - 2Awi -ia + z'/Jy, (40) 

C c = -2<tai +if3uj 1 + 2Alu 1 + ia-i/3y. (41) 

Both A and B are therefore functions of the amplitude and detuning of the excitation 
drive (F and 6, in the rotating frame) and the axialization drive (e and 2A), the laser 
parameters (a and (3) and the trap frequencies (to c and toi). Note that the ratio A/B 
is independent of the driving force amplitude, F, and remains the same as F tends 
to zero. As expected we find that in the absence of axialization, when e = 0, the 
amplitude of one of the components of the motion, B, is also zero. 

A and B are complex quantities and the amplitude and phase response of the 
ion's motion relative to the excitation drive is obtained by taking their moduli and 
arguments respectively. Figure 7 shows a typical contour plot of the amplitude 
response \A\ as a function of axialization and excitation drive detuning (in the 
laboratory frame) for a very small axialization amplitude. The plotted parameters 
are identical to those used in figure 5(a). This plot corresponds to the response of the 
motion being driven, and hence is the response that can be observed experimentally. 
It is clearly evident from the plot that \A\ is only large when the excitation drive 
detuning is close to zero, and there is no indication of a second resonance frequency. 
An equivalent plot of |B| would correspond to the response of the component of the 
motion that is not being directly driven (i.e. the counter-rotating component) and 
would be difficult to measure experimentally When the axialization amplitude is very 
small, \B\ is found to be almost zero regardless of the axialization and excitation 
drive dctunings. Essentially this confirms that if there is no coupling between the 
two motions then driving one of them will have no effect upon the other. These 
results validate the earlier assertion that the components of the motion with the extra 
frequencies observed in the model have very low amplitudes when the axialization 
drive amplitude is small or its detuning is large. 

Also shown in figure 7 is Arg(A), the motional phase of the driven component 
of the motion relative to the excitation drive. The parameters used for this plot 
correspond to the dotted line marked on the contour plot of \A\. There is a transition 
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from being in phase with the excitation drive to lagging tt radians behind it as the 
drive is tuned through resonance. This is identical to the phase response of a single 
classical oscillator relative to a driving force and is in agreement with measurements 
made in the absence of axialization [38]. 

Figure 7 shows the amplitude response of the ion's motion when the axialization 
drive amplitude is large. The parameters used to generate these plots are the same 
as those used in figure 5(d). If the axialization drive is nearly resonant then the 
motion can be driven at two distinct frequencies. The amplitude response |A| to the 
excitation drive at each of these frequencies is similar — as expected for a strongly 
coupled system. Furthermore, driving one component of the motion, A, leads to a 
build up of amplitude in the other component, B, when close to resonance. This is 
a direct indication of the efficacy of the axialization technique, since a damping force 
on one component of the motion will likewise lead to a reduction in amplitude of the 
other component. 

The phases of the two components of the ion's motion relative to the excitation 
drive, given by Arg(A) and Arg(B), are also plotted in figure 7. Once again the 
parameters used are indicated by the dotted lines on the corresponding amplitude 
contour plots. In addition to the 7r phase shift observed in the case of weak axialization, 
Arg(A) now exhibits a second phase shift at the position of the second resonance 
frequency. The size of this second phase shift is found to tend towards it as the 
axialization amplitude is increased. 

Arg(.B) also undergoes two phase shifts. Unlike the component of the motion 
A, which always lags behind the excitation drive and so always has a relative phase 
between — tt and zero, the relative phase of B can take any value. 

The predicted phase behaviour of the measurable component of ion motion, 
Arg(A), is in qualitative agreement with recent experimental results. 

6. Discussion 

We have presented an analytical model for the axialization of laser-cooled ions in 
a Penning trap. This model allows the calculation of the frequencies of motion of 
the ions and the behaviour is found to be different in each of three distinct regimes 
corresponding to weak, intermediate and strong coupling due to axialization. The 
first and last of these regimes are easily accessed experimentally and the results of 
the model are in good agreement with measured frequencies. The model can also be 
used to study the amplitude and phase response of the axialized ions relative to an 
oscillating dipolar electric potential, yielding results that are in qualitative agreement 
with measured data. 

The model demonstrates that axialization can be used as a tool to enhance the 
efficiency of laser cooling of ions trapped in Penning traps and allows the calculation of 
the range of parameters for which such an improvement may be obtained. By helping 
to overcome this significant experimental difficulty, axialization can make the Penning 
trap a viable alternative to the RF trap for quantum information processing and other 
experiments where precise control of a single trapped atomic particle is required. 
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Figure 7. Left: Contour plots showing the amplitude response of a calcium ion to 
an external dipolar drive near the magnetron frequency when (e/aJi) 2 = 0.01A/ 2 
and (e/oji) 2 = 100M 2 , corresponding to cases (a) and (d) in figure 5. \A\ is 
the response of the component of the motion being directly driven, whilst \B\ 
is the response of the indirectly driven component. Note that the plot of \B\ is 
not shown for the case when (e/o>i) 2 = 0.01M 2 , since it is very close to zero for 
all values. The scale is the same for each plot, but the absolute amplitudes are 
dependent on the driving force and are therefore arbitrary. Right: Plots of the 
phase of each component of the ion motion relative to the excitation drive. The 
parameters used correspond to the dotted line on the contour plots. In each case 
u) c = 380x27rkHz, u) 1 = 165x2tt kHz, a/j3 = 100x2tt kHz, and |M| 2 = 0.01 kHz. 
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Appendix A. Solution of the equations for laser cooling in the presence of 
axialization 

Substituting trial solution (21) into (20) yields 



Aui 2 - A(ip + <j c - 2u) r )ui 

2 



+e" 



+ A ( -1 +ia + uj 2 - uj t (uj c + i/3) ) - eB* 
Buj* 2 + B(if3 + uj c - 2ui r )uj* 



B[^-+ia + Lu 2 - uo r {uo c + i(3) 



eA* 



= 0. 



(A.l) 



It is straightforward to show that if c • f{x) + c* ■ g(x) = 0, where c has non-zero real 
and imaginary parts, then f{x) and g(x) must both be equal to zero. Therefore 



u> — (i(3 + uj c — 2uj r ) uj 

2 



B 



+ I + ia + w r — w r (w c + i(3) 



uj* 2 + (if3 + uj c - 2Lu r ) uj* 

( U) 2 

+ I + ia + uj 2 - Lu r (to c + i/3) 



eB* 



eA* = 0. 



(A.2) 



(A.3) 



It was shown in section 4 that in order to be in a frame in which the axializing 
quadrupolar drive is static it is necessary to have uj r = ui a /2. A is therefore defined 
to be half the detuning of the axialization drive from the cyclotron frequency, so that 
uj a = w c + 2A. Thus the frame rotation frequency is w r = lu c /2 + A. In addition to this 
the variable S is defined to be the (complex) shift in the frequencies from u)\ due to 
the presence of laser cooling, the axializing drive and the change to the frame rotation 
frequency uj r , so that u) = u>\ + 5. The equations for ui can therefore be rewritten as 



A 



B 



28u!i — i(3u)\ + 2Acji 
. 2a - Qui, 



+ (A + 5) - i/3 (A + 5) 



25*lui + ifiuii — 2Awi 
2a - pu>, 



+ (A - S*) - ip (A - S* 



eB* = 



- eA* = 0. 



(A.4) 



(A.5) 



The (A ± 6) and P(A±S) terms are small relative to those containing ui\ . Ignoring 
these, and taking the complex conjugate of the latter equation, the pair of equations 
becomes 

. 2a — Pui c 



(2<Jwi - iPu)\) + 2Acji + i 



eB* 







(A.6) 
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B* 



(25o;i - ipui) - ( 2Auji + i 



2a — j3u! c 



Solving (A. 7) for A and substituting into (A. 6) gives 



B* 



(2fci - i(3uJif - { 2Acji + i 



2a — [3lo c 



e*A = 0. 



= 0. 



(A.7) 



(A- 



This can be simplified by defining M = (2a — (3oj c ) /2ui, a term related to the laser 
cooling parameters a and j3. In order to find the real and imaginary parts of this 
expression the substitution 5 = Sq + 170 is also made, so that 

= 0. (A.9) 



(25 + 2i l0 - i/3) 2 - (2A + iMf 



All of the variables arc now real. Evaluating the squared terms and separating out 
real and imaginary parts gives 

I |2 



A8l - 4 7 2 ~p 2 + 4/3 7o - 4A 2 + M 







and 



8(5 7o ~ 4<5 /3 - 4AM = 0. 
The latter equation gives the cooling rates 
(i AM 

70= 2 + W 
Substituting 70 into (A. 10) then yields 



51 



M 2 



4 W 2 



A 2 M 2 

-T3-=°. 



45, 



(A.10) 
(A.ll) 
(A.12) 

(A.13) 







which, if N is defined to be equal to the term in brackets, leads to the quartic equation 



St N5£ 
This has four roots 
So 



A 2 M 2 



1 



±^^N±^/N 2 + A 2 AP. 
V2 



(A.14) 



(A.15) 



Since <5o is by definition real, the positive sign must be taken for the inner square root 
and so two of the solutions can be eliminated. Indeed these additional complex roots 
are in fact duplicated solutions. If we use them to determine 7 o we obtain identical 
expressions for S = 5q + i"fo- 



Appendix B. Calculation of motional amplitudes for a driven, laser-cooled 
ion in the presence of axialization 

A dipolar field oscillating at frequency u>d in the laboratory frame can be thought of 
as a pair of static dipolar fields rotating in opposite senses at the oscillation frequency, 
since 

2Fcosu; d t = F (e iuJdt +e~ luJdt ) . (B.l) 
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If the driving frequency is close to the magnetron frequency, so that lo^ = uj m + 5' , 
then in a frame rotating at a frequency uj i = (cj c /2) + A the driving force is given by 



F 



(B.2) 



The first term here is nearly resonant with the frequency of the magnetron motion 
in this frame (the detuning in the rotating frame being 5' — A, henceforth defined to 
be S). The second term will in general be far from resonance and can be ignored. 
uj is defined to be the frequency of the driving field in the rotating frame, so that 
uj = uj c \ — (w c /2) — A = —oj\ + 5. 

Adding in this driving force into the equation of motion in the rotating frame for 
an ion being laser cooled with axialization (20) leads to 

v + (/? — iuj c + 2iuj r )v — ( -2- + ia + u% — u) r (oj c + i(3) ) v + ev* 



= Fc lut . (B.3) 

Steady state solutions corresponding to the situation once the system has reached 
equilibrium will have the form of elliptical motion with a frequency lu, such that 



Be 



-IUJ* t 



(B.4) 



Here to is real and is the known applied frequency rather than a variable. Proceeding 
in an identical manner to before leads to the equivalent of (A. 6) and (A. 7): 



.4 



B* 



— 2&ji + ifiuoi — 2Awi — ia + 

— 2Su!i + i(3uJi + 2At^! + ia — 



eB* = F, 
e*A = 0. 



(B.5) 



(B.6) 



The first of these equations is for the component rotating in the same sense as the 
driving force (in this case the direction of the magnetron motion), with the second 
equation describing the component of rotation in the opposite sense. Bearing in mind 
that 5 is now a defined constant detuning, the constant terms in square brackets can 
be grouped together as C m and C c respectively. This gives simply 

AC m + eB* = F, 



B*C c +e*A = 0. 
These equations are easily solved for A and B, yielding 
A Fa 



B = 



C m C c - \e[ 
Fe 



(B.7) 
(B.8) 

(B.9) 
(B.10) 
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